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1.  SUMMARY 


1.1.  OBJECTIVE  AND  ACCOMPLISHMENTS 

The  overall  objective  of  the  proposed  research  was  to  (a)  fundamentally  and 
quantitatively  understand,  based  on  first  physical  principles,  the  micromechanics  of 
penetration  into  and  damage  accumulation  in  cemented  peuticulate  targets  subject  to 
multiple  impact:  and  (b)  translate  this  understanding  into  a  predictive  code  by  using 
contact  and  damage  laws  thus  derived  in  a  DEM  routine.  The  technical  approach  was 
to  integrate  experiments,  analyticed  micromechanics,  and  numerical  modeling  into  a 
concerted  effort.  The  project  was  to  be  conducted  jointly  with  the  University  of  Rhode 
Island.  The  Petrophysical  Consulting  part  of  the  project  was  analytical  micromechanics. 

To  achieve  this  objective  we  continued  our  previous  effort  (AFOSR  Contract  F49620- 
96-1-0394)  and  provided  an  elastic  description  of  completely  cemented  (no  void)  or 
partially  cemented  granular  material  where  the  elastic  properties  of  the  grains  and 
cement  differ  from  each  other.  This  new  result  provides  a  solution  to  the  problem  of 
elastic  composite  with  contacting  inclusions  —  the  problem  that  has  not  been  solved 
before.  This  solution  gives  a  basis  for  investigating  the  plastic  deformation  and  failure  of 
cemented  composite  targets  such  as  concrete  and  granite. 

Our  second  achievement  was  in  analyzing  natural  rock  materials  using  acoustic 
pulses  and  deriving  the  strength  of  such  materials  from  their  elastic  moduli.  These 
results  are  relevant  to  describing  the  strength  and  damage  resistance  of  such  target 
geomaterials  as  concrete  and  granite. 

The  main  relevanee  to  the  Air  Force  mission  is  through  a  quantitative  description  of 
elastic,  non-elastic,  and  damage  accumulation  properties  of  cemented  geomateiials 
subject  to  penetrator  impact. 


1.2.  PERSONNEL  SUPPORTED 

Senior  Personnel:  Dvorkin,  Jack,  PI  (US);  Geuy  Mavko  (US);  eind  Amos  Nur  (US). 

1.3.  PUBLICATIONS 

1.  Jack  Dvorkin,  Jim  Berryman,  Amos  Nur,  and  Mickaele  Le  Ravalec,  1998,  Effective 
moduli  of  particulates  with  elastic  cement.  Mechanics  of  Materials,  in  press. 

2.  Jaek  Dvorkin,  and  Ivar  Brevik,  1998,  Diagnosing  high-porosity  seindstones:  Strength 
and  permeability  from  porosity  and  velocity.  Geophysics,  in  press. 

3.  Per  Avseth,  Jack  Dvorkin,  Gary  Mavko,  and  Johannes  Rykkje,  1998,  Diagnosing 
high-porosity  sands  for  reservoir  charaeterlzatlon  using  sonic  and  seismic,  1998  SEG 
Annual  Conference,  New  Orleans,  RP  2.5,  1024-1027. 

1 .4.  INTERACTIONS  AND  TRANSITIONS 

•  1998  Stanford  Rock  Physics  Meeting;  Stanford,  CA,  1998:  Natural  cemented  roeks. 

•  1998  SEG  Annual  Conference,  New  Orleans:  Diagnosing  high-porosity  sands  for 
reservoir  characterization  using  sonic  and  seismic. 

•  ASCE  1998  Meeting,  La  Jolla:  Effective  Moduli  of  Particulates  with  Elastic  Cement. 

1.5.  PATENTS  AND  INVENTIONS 

No  patents  or  inventions  resulted  form  this  effort  (see  DD  882  attached). 
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2.  TECHNICAL  PART 


In  this  part  the  technical  results  are  presented  In  the  form  of  journal  articles.  The 
numeration  of  formulas  and  figures  Is  Internal  for  every  section. 

2. 1.  EFFECTIVE  MODULI  OF  PARHCULATES  WITH  ELASTIC  CEMENT 
In  collaboration  with  Jim  Berryman,  Amos  Nur,  and  Mlckaele  Le  Ravalec 


ABSTRACT 

We  give  an  approximate  method  for  calculating  the  effective  elastic  moduli  of  a  close 
random  pack  of  Identical  elastic  spheres  whose  pore  space  Is  partially  or  eompletely  filled 
with  elastic  cement.  To  construct  the  solution  we  start  with  the  uncemented  pack  and 
add  small  amounts  of  cement  uniformly  at  every  grain  eontact.  The  effective  moduli  of 
such  an  aggregate  depend  on  those  of  the  grains  and  cement  (which  are.  generally, 
different),  and  on  the  amount  of  cement.  Moduli  for  grain/cement  mixtures  are 
available  from  the  contact  cementation  theory.  Next  we  introduce  an  isotropic  elastic 
body  with  void  inclusions  (e.g.,  spherical)  whose  porosity  Is  the  same  as  the  porosity  of 
the  pack  with  contact  cement.  The  moduli  of  the  matrix  are  as  yet  unknown.  We  find 
them  by  assuming  that  the  effective  moduli  of  the  elastic  body  with  inclusions  are  equal 
to  those  of  the  pack  with  contact  cement.  Then  an  appropriate  effective  medium  theory 
(e.g.,  self-consistent  or  differential  effective  medium  approximation)  provides  two 
equations  for  the  bulk  and  shear  moduli  of  the  matrix.  Finally,  we  use  the  chosen 
effective  medium  theory  to  calculate  the  moduli  of  the  elastic  body  with  the  matrix  thus 
defined  whose  inclusions  are  now  filled  with  the  cement.  We  assume  that  these  are  also 
the  moduli  of  the  pack  whose  pore  space  is  completely  filled  with  the  cement.  To 
caleulate  the  elastic  moduli  of  the  pack  whose  pore  space  is  partially  filled  with  the 
cement,  we  assume  that  they  are  those  of  the  elastic  body  with  the  matrix  thus  defined 
and  with  two  types  of  inclusions  -  void  and  filled  with  cement  (the  porosity  of  this  body 
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and  of  the  cemented  pack  are  the  same).  This  result  gives  one  solution  to  a  general  smd 
previously  unresolved  problem  of  contacting  elastic  inclusions  embedded  in  an  elastic 
matrix.  It  provides  a  transition  from  high-porosity  cemented  to  completely  cement-filled 
granular  materials.  The  solution  accurately  predicts  the  compressional  elastic  modulus 
(the  product  of  the  bulk  density  and  P-wave  velocity  squared)  from  wave-propagation 
experiments  on  epoxy-cemented  glass  beads,  ice-cemented  Ottawa  sand,  and  in-situ 
acoustic  data  in  natural  sands  cemented  by  gas  hydrate. 

INTRODUCTION  AND  PROBLEM  FORMULATION 

When  estimating  the  effective  stiffness  of  many  composites,  concrete  and  natural 
rock  among  them,  one  often  faces  the  problem  of  calculating  the  moduli  of  an  elastic 
matrix  with  elastic  inclusions.  This  problem  has  been  extensively  addressed  in  the  past. 
The  existing  solutions  are  analyzed  and  summarized  by,  e.g.,  Christensen  (1991), 
Zimmermein  (1991),  Wang  and  Nur  (1992),  Nemat-Nasser  and  Hori  (1993),  and  Benyman 
(1995). 

Accurate  solutions  are  available  for  small  inclusion  concentrations  and  various 
inclusion  shapes.  To  obtain  an  estimate  for  high  concentrations,  one  may  find  the  self- 
consistent  (SC)  or  differential  effective  medium  (DEM)  approximations  to  be  useful  and 
reasonably  accurate.  Some  fundamental  solutions  for  composites  with  concentrated 
interacting  inclusions  belong  to  Frankel  and  Acrivos  (1967),  Goddard  (1977),  and  Chen 
and  Acrivos  (1978). 

Ju  and  Chen  (1994)  show  how  to  calculate  the  effective  elastic  moduli  of  two-phase 
composites  containing  randomly  dispersed  sphericed  inhomogeneities.  However,  their 
solution  does  not  treat  the  importemt  case  of  spheres  in  direct  contact  and,  therefore, 
fails  to  predict  experimentally  measured  elastic  moduli  in  such  composite  (see  section 
Examples  and  Appendix  2). 
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An  alternative  to  direct  estimates  is  the  bounding  of  the  overall  moduli  of  a 
composite  (Hashin  and  Shtrikman,  1963).  Nemat-Nasser  et  al.  (1993)  used  generalized 
Hashln-Shtrikman  variational  principles  to  provide  accurate  bounds  for  composites 
with  periodic  microstructure.  They  also  introduced  a  "modified  equivalent  inclusion 
method"  for  directly  estimating  the  moduli  of  such  composites.  Tlie  results  fall  between 
the  upper  and  lower  bounds  and  are  supported  by  experimental  data  on  the  shear 
modulus  of  an  incompressible  matrix  reinforced  by  rigid  spherical  particles. 

In  composites  with  concentrated  elastic  inclusions,  direct  interaction  between  the 
inclusions  through  the  thin  layers  of  the  matrix  is  very  important  (Goddard,  1977).  The 
contact  cementation  theory  (CCT)  of  Dvorkin  et  al.  (e.g.,  1994)  highlights  this 
conclusion.  It  allows  one  to  calculate  the  effective  elastic  properties  of  a  particulate 
aggregate  where  cement  is  located  eiround  the  grain  contacts  (Figure  la). 

CCT  predicts  that  even  small  volumes  of  contact  cement  (which  may  be  soft  as 
compared  to  the  grains)  act  to  strengthen  the  aggregate's  elastic  moduli,  thus  turning 
the  nonlinear  problem  of  contact  analysis  into  a  simpler  lineeir  one  such  as  that  found 
typically  in  composites  analysis.  The  initial  filling  of  the  contact  gaps  is  most  important 
—  the  effect  of  any  additional  cement  placed  around  this  initial  filling  is  relatively  small. 
This  result  is  supported  by  experiments  (Figure  2).  Experiments  also  show  that  even  by 
filling  the  entire  pore  space  of  a  sphere  pack  one  cannot  achieve  the  drastic  relative 
stiffness  increase  for  which  the  small  volumes  of  the  contact  cement  are  responsible. 

The  assumptions  used  in  CCT  make  its  use  appropriate  only  for  high-porosity 
particulates  where  cement  is  located  in  the  immediate  vicinity  of  the  contacts.  As  is, 
CCT  cannot  be  used  to  estimate  the  elastic  constants  of  a  zero-  or  low-porosity  aggregate 
where  cement  fills  the  whole  pore  space  (Figure  lb)  or  large  portions  of  it. 
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Therefore,  a  question  naturally  arises:  How  do  we  calculate  the  effective  elastic 
moduli  of  particulates  where  cementation  extends  far  beyond  immediate  graln-to-graln 
contacts?  Consider,  for  example,  a  close  random  pack  of  identical  glass  beads  where  the 
pore  space  is  completely  filled  with  epoxy.  This  aggregate  is  a  case  of  contacting 
spherical  inclusions  (glass)  embedded  in  the  epojqr  matrix.  From  CCT  we  know  that  the 
filling  (epoxy)  located  at  the  contacts  contributes  most  to  the  stiffness  of  the  material. 
Then  it  is  clear  that  a  direct  use  of  an  Inclusion  theory  (such  as  SC  or  DEM)  is 
inappropriate  here  because  it  does  not  account  for  the  special  properties  of  cemented 
contacts.  Contact  cement  and  pore-filling  cement  (Figure  Ic)  contribute  differently  to 
the  stiffness  of  the  material. 


a 


b 


ffl 


Figure  1.  a.  Grains  with  contact  cement,  b.  Grains  with  pore  space  completely  filled,  c. 
Contact  cement  and  pore-filling  cement. 


To  treat  this  problem  our  approach  is  to  obtain  a  solution  to  the  concentrated 
inclusion  problem  (the  bulk  and  shear  moduli  of  a  close  random  pack  of  identical  elastic 
spheres  whose  pore  space  is  partially  or  completely  filled  with  elastic  cement)  where  the 
(previously  ignored)  special  properties  of  the  contact-cement  part  of  the  matrix  are 
recognized  and  included.  The  problem  posed  is  Important,  not  only  methodologically, 
because  many  cemented  geomaterials,  such  as,  e.g.,  concrete,  have  inclusion 
concentrations  close  to  the  upper  limit,  small  void  concentration  (porosity)  and  large 
amounts  of  pore-filling  cement. 
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Cement  Saturation  of  Pore  Space 


Figure  2.  CCT  theoretical  predictions  (curves)  and  experimental  data  (symbols)  for  the 
compressional  elastic  modulus  from  elastic  wave  propagation  experiments.  The  cemented 
samples  were  prepared  by  (a)  adding  epo^^  to  a  close  random  pack  of  identical  glass  beads 
(Yin,  1993),  and  (b)  adding  water  to  Ottawa  sand  and  then  freezing  it  (Tutuncu,  1997). 
Both  epoxy  and  water  are  wetting  fluids  and,  when  they  solidify,  act  as  contact  cement  in 
the  form  of  pendular  rings.  The  modulus  is  plotted  versus  the  cement  saturation  of  the 
pore  space  of  the  uncemented  glass  bead  and  sand  packs. 


SOLUTION 

To  solve  this  problem,  we  assume  that  at  some  high  porosity  where  CCT  is  still 
applicable  (Figure  3,  point  B),  the  effective  bulk  and  shear  moduli  of  the  cemented 
material  are  those  of  a  homogeneous  isotropic  matrix  with  empty  spherical  inclusions  of 
the  same  porosity.  Next  we  use  am  effective  medium  theory,  for  example,  SC  (Berryman, 
1980)  or  DEM  (Norris,  1985),  to  express  these  effeetive  bulk  and  shear  moduli  through 
the,  as  yet  unknown,  elastic  moduli  of  the  h5^othetical  matrix  material.  Then  the 
elastic  moduli  of  the  matrix  material  can  be  found  from  the  two  resulting  equations. 

To  find  the  elastic  moduli  of  the  sphere  pack  where  the  entire  pore  space  is  filled  with 
cement,  we  assume  that  they  are  the  same  as  those  of  the  h5q)othetical  matrix  whose 
inclusions  are  now  completely  filled  with  the  cement  (Figure  3,  point  C).  Those  moduli 
can  be  computed  using  the  chosen  effective  medium  theory. 

Finally,  we  find  the  moduli  of  the  cemented  sphere  pack  of  a  non-zero  porosity  by 
assuming  that  they  are  those  of  the  cement-filled  hypothetical  matrix  into  which  empty 
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spherical  inclusions  are  embedded  (Figure  3,  point  D).  The  porosities  of  this  new  body 
with  inclusions  and  of  the  cemented  pack  are  the  same. 

The  mathematical  expressions  for  this  cdgorlthm  are  given  in  Appendix  1. 

Of  course,  the  inclusions  do  not  have  to  be  spherical  (results  for  ellipsoidal-shaped 
inclusions  are  easily  obtained  within  the  same  approach),  cind  effective  medium  theories 
other  than  SC  emd  DEM  could  be  used. 


Figure  3.  Proposed  solution  scheme.  "A",  uncemented  sphere  pack  of  zero  stiffness.  "B", 
pack  with  small  volume  of  contact  cement  whose  moduli  are  calculated  from  CCT;  the 
moduli  of  the  hypothetical  matrix  material  are  found  from  SC  or  DEM.  "C",  completely 
cemented  matrix  and  sphere  pack.  "D",  matrix  and  sphere  pack  of  non-zero  porosity.  The 
graph  (compressional  modulus  versus  porosity)  is  given  for  epoxy-cemented  glass  beads. 


EXAMPLES 

In  the  examples  below,  the  compressional  modulus  M  is  caleulated  from  the 
measured  P-wave  velocity  and  bulk  density  p  as  M  =  pVp.  In  the  laboratory 
experiments  presented  here,  the  shear-wave  velocity  has  been  recorded  as  well.  We  do 
not  use  these  experimental  results  because  the  quality  of  the  shear-wave  data  was 
(traditionally)  poor  due  to  the  interference  of  converted  P  -waves.  The  wave-length  in  all 
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the  experiments  is  much  larger  that  the  individual  grain  size,  which  justifies  the  effective 
medium  approach  to  the  problem. 

The  elastic  properties  of  materials  (glass,  quartz,  epoxy,  and  ice)  used  in  the 
calculations  are  given  in  Table  1 . 


Table  1.  Elastic  moduli  of  materials  used  in  calculations. 


Materiail 

Shear  Modulus  (GPa) 

Bulk  Modulus  (GPa) 

Glass 

30.34 

45.55 

Epo;^^ 

2 

6.8 

Quartz 

45 

36.6 

Ice 

3.53 

8.57 

Gas  Hydrate 

2.4 

5.6 

LABORATORY  EXPERIMENTS.  The  elastic-wave  velocities  have  been  measured  in 
glass  bead  packs  cemented  with  epo^qr  by  Yin  (1993),  and  in  water-saturated  and  frozen 
pure-quartz  Ottawa  sand  by  Tutuncu  (1997).  Both  liquid  epojqr  eind  water  are  wetting 
fluids  and  tend  to  accumulate  at  grain  contacts.  Therefore,  at  high  porosity  CCT  is 
applicable  to  estimating  the  effective  elastic  moduli  of  the  aggregates  under  examination. 

At  low  porosity  we  use  the  above-described  solution  with  SC  (Berryman,  1980)  and 
DEM  (Norris,  1985)  approximations  for  spherical  inclusions.  In  both  cases  our 
theoretical  estimates  match  the  experimental  data  well  (Figure  4). 

An  adjustable  parameter  in  the  solution  is  porosity  above  which  CCT  is  used  and 
below  which  the  new  solution  is  employed.  Our  solution  is  stable  with  respect  to  this 
parameter  (see  Figure  5  for  epo^-cemented  glass  beads)  and,  therefore,  this  transition 
porosity  can  be  arbitrarily  chosen  within  the  reasonable  limits  where  CCT  is  applicable. 

It  would  be  more  satisfying  if  all  of  the  new  solution  curves  Intersected  the  zero 
porosity  axis  at  the  same  modulus  value,  but  instead  we  find  that  the  present 
formulation  of  the  theory  produces  a  set  of  curves  that  are  shifted  (more  or  less)  in 
parallel  depending  on  their  point  of  intersection  with  the  CCT  curve.  This  result 
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suggests  that  the  current  formulation  might  be  improved  using  a  more  sophisticated 
effective  medium  method. 


Figure  4.  Compressional  modulus  versus  porosity  for  (a)  epoxy-cemented  glass  beads,  and 
(b)  ice-cemented  Ottawa  sand.  Filled  symbols  are  from  the  experiments.  Curves  are  from 
theory.  Open  S5mibols  show  the  results  of  applying  effective-medium  theories  to  the  grain- 
cement  mixture  at  zero  void  concentration. 


The  most  important  result  of  our  new  solution  is  the  elastic  moduli  of  cemented 
spheres  at  zero  void  concentration.  In  order  to  show  that  it  is  crucial  to  account  for  the 
specific  geometry  (contacting  spheres)  of  the  aggregate,  we  compare  this  solution  to 
various  two-component  effective  medium  results.  We  use  the  Hashin-Shtrikman  (1963) 
bounds  (Upper  HS  and  Lower  HS),  SC  (Berryman,  1980),  Lower  DEM  (Norris,  1985),  and 
the  Ju  and  Chen  (1994)  solution  (JC)  to  directly  mix  cement  and  spheres  at  zero  void 
porosity.  The  results  are  shown  in  Figure  4  and  given  in  Table  2.  Only  the  Lower  HS 
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values  appear  to  well  match  the  experimental  data  and  our  new  solution.  Tliis  is 
consistent  with  the  geometry  of  the  aggregates  under  examination  where  softer  cement 
envelopes  stiffer  grains.  All  other  effective  medium  theories  fail  to  match  the 
experimental  data. 


Figure  5.  Stability  of  the  new  solution.  Compressional  modulus  of  epoxy-cemented  glass 
beads  versus  porosity.  Bold  curve  are  from  CCT,  thin  lines  are  from  our  new  solution  (SC). 
These  lines  start  at  a  varying  transitional  porosity. 


Table  2.  Appl3dng  various  effective-medium  theories  at  zero  porosity. 


Compressional  Modulus  (GPa) 

Glass-Epoxy 

Quartz-lee 

Experiment 

24.2 

34.0 

This  Paper 

25.4 

35.7 

SC 

39.0 

44.0 

Lower  DEM 

31. 

40.0 

Upper  HS 

46.3 

51.9 

Lower  HS 

26.0 

32.7 

JC 

43.1 

42.7 

IN-SITU  MEASUREMENTS.  Gas  hydrate  is  an  ice-like  substance  resulting  from  a 
reaction  between  methane  and  water  in  a  favorable  temperature  and  pressure  range. 
Hydrates  discovered  in  nature  often  fill  the  pore  space  in  natural  rocks  and  cement  the 
rock  grains.  An  example  presented  here  is  based  on  in-situ  velocity,  porosity,  and  gas- 
hydrate  saturation  measurements  conducted  in  a  well  drilled  through  sandy  intervals  in 
the  North  Slope  of  Alaska  (Mathews,  1986). 
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In  Figure  6  we  give  the  porosity  of  the  rock  frame  (without  gas  hydrate),  volumetric 
fraction  of  gas  hydrate  in  the  pore  space  (from  electrical  measurements),  and  the 
compressional-wave  velocity  as  measured  in  situ  by  well  logging  tools. 

The  acoustic  data  Is  compared  to  our  theoretical  prediction  (Figure  6c).  In  this  case 
we  used  our  new  solution  with  SC  approximation.  The  match  between  the  data  and 
theory  is  very  good. 


Figure  6.  Porosity  (a),  gas  hydrate  concentration  (b),  and  compressional-wave  velocity  (c) 
versus  depth  in  a  gas  hydrate  well. 


CONCLUSION 

We  calculate  the  effective  moduli  of  an  elastic  matrix  with  elastic  inclusions  at  high 
concentration  by  modeling  it  as  a  cemented  particulate  aggregate.  This  approach  allows 
us  to  account  for  the  important  role  of  the  contact  cement  part  of  the  matrix  and 
address  the  problem  of  directly  contacting  inclusions.  The  latter  problem  cannot  be 
approached  using  traditional  solutions  where  the  inclusions  are  considered  absolutely 
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rigid  (e.g.,  Frankel  and  Acrivos,  1967)  and  separated  by  a  finite-thickness  layer  of  the 
matrix. 

The  solution  is  a  combination  of  the  contact  cementation  theory  emd  an  effective 
medium  theory  such  as  the  self-consistent  or  differential  effective  medium 
approximation.  The  solution  is  approximate  because  it  uses  a  number  of  assumptions. 
Nevertheless,  its  predictions  are  confirmed  by  experimental  data  on  cemented  glass 
beads,  sand,  and  natural  rocks.  We  are  able  to  accurately  model  the  dynamic 
compressional  moduli  of  these  aggregates.  None  of  the  existing  effective-medium 
solutions  can  adequately  match  the  experimental  results. 

One  legitimate  criticism  of  the  present  formulation  of  the  theory  is  that  neither  of  the 
two  effective  medium  theories  used,  although  well-established  and  relatively  easy  to  use, 
permits  the  user  to  control  the  microstructure  of  the  resulting  composite  in  any 
substantial  way.  In  particular,  both  approaches  when  applied  in  Steps  3  and  4  (see 
Appendix)  will  surely  violate  our  assumed  microstructure  for  intermediate  porosities  by 
effectively  placing  voids  in  the  grains  (where  we  do  not  want  them)  as  well  as  in  the 
cement  material  (where  we  really  do  wemt  them  to  be). 

This  effect  will  be  comparatively  small  when  the  cementing  material  has  similar  or 
the  same  elastic  parameters  as  the  grains  as  is  the  case  for  our  glass-glass  and  quartz- 
quartz  examples.  But  it  may  have  a  substantial  effect  on  the  predicted  properties  in  the 
cases  where  the  ratio  of  grain  to  cement  properties  is  large  as  is  the  case  for  glass-epo^ 
and  quartz-ice.  The  expected  bias  that  would  result  from  the  undesired  voids  in  the 
grains  is  that  the  predicted  effective  values  of  the  gredn/cement  mixture  at  zero  porosity 
must  contain  a  higher  proportion  of  grain  material  in  the  total  volume  than  that 
present  in  the  CCT  calculations.  The  result  will  be  a  bias  towards  higher  them  realistic 
values  at  zero  porosity  and  this  is  exactly  what  we  observe  for  the  SC  theory  in  the 
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glass-epoxy  and  quartz-ice  examples.  From  this  point  of  view,  the  present  results 
should  be  considered  preliminary;  we  need  to  continue  searching  for  relatively  simple 
effective  medium  theories  that  will  permit  us  to  gain  control  over  the  microstructure 
being  modeled. 

Still,  this  solution  cem  be  directly  applied  to  estimating  the  stiffness  of  particulates 
such  as  concrete  where  the  grains  touch  each  other  emd  the  void  fraction  is  small. 
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APPENDIX  1:  SOLUTION  ALGORITHM 


Step  1.  Contact  Cementation  Theory  fCCT).  The  effective  compressional  (A/^ct),  bulk 
(Kccr),  and  shear  moduli  of  a  random  dense  pack  of  identical  spheres  cemented 

at  their  contacts  are  (Dvorkin  and  Nur,  1996); 


where  £ind  are  the  compressional  and  shear  moduli  of  the  cement  respectively; 
^0  0. 36  is  the  porosity  of  the  uncemented  pack;  eind  n~S.5  is  the  average  number  of 

contacts  per  grain. 

Parameters  and  are: 


S„  =  A„(A„)a'  +  5„(A„)a  +  C„(A„),  A.(A„)  =  -0.024153- A„-'•'^^ 

=  0.20405 •  C„(AJ  =  0.00024649  •  A„-* 

=  A,(A^,  v)a^  +  fi,(A^,  v)a  +  C,(A„  v), 

A,(A„  V)  =  -10''  •  (2.26 v'  +  2.07V  +  2.3) • 

V)  =  (0.0573v^  +0.0937  v+ 0.202) -A,®®'’""' 

C,  ( A  „  V)  =  1 0-^  •  (9. 654  +  4. 945  V  +  3. 1)  •  A  +o.4oi  t  v-i.8i86 . 

_2G,(l-v)(l-vJ 
"  ttG  1-2v,  ’  "  tcG’ 

where  G  and  V  are  the  shear  modulus  and  the  Poisson's  ratio  of  the  grain  material, 
respectively:  and  is  the  Poisson's  ratio  of  the  cement.  Parameter  OC  is 


a  =  2[ 


^0  ~  ^  -jO.25. 

3n(l-0o)  ’ 


(j)  is  the  porosity  of  the  cemented  aggregate  which  is  smaller  than  because  of  cement 
deposited  in  the  pore  space  of  the  sphere  pack. 

Step  2.  Self-Consistent  Approximation  fSCL  If  the  moduli  and  G^  of  the 
hypothetical  matrix  materieil  were  known,  then  the  self-consistent  approximation  for 
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spherical  inclusions  of  void  with  concentration  (j)  predicts  that  ^ccr  given 

by  the  coupled  equations  (Berryman,  1980): 


1-^  ,  0 


^CCT  ^  ^CCT  ^5  ^  ^CCT  ^  ^CCT 

1  _  ^  ~  ^  0  ^  _  ^ccr(^^ccT  ^^ccr) 

^ccr  Z  +  Z  Z  +  2G(;.(^) 

These  equations  have  to  be  solved  iteratively  to  find  and  G^cr  when  and 
Gj  are  known.  But,  since  K(^(^  and  Gccr  already  given  in  Step  1,  we  may  solve  the 
two  equations  Instead  for  and  G^ .  Solving  in  this  direction,  the  equations  decouple 
and  provide  explicit,  unique  results  for  and  G, . 

Step  2'.  Differential  Effective  Medium  Approximation  (DEMI.  In  this  case  euid  G^ 
can  be  found  fi*om  the  following  two  equations  (Norris,  1985): 

^CCT  _  /I  r ^  ^  ■^i^CCT  !  ^CCr  _  ^CCT  _ 0.75  +  _ 


where 

T  3(1- 5  vj  13^.-2G, 

^  4(1+ V,)’  ^  2  3K^+G/ 


Norris'  result  can  be  rewritten  to  show  explicitly  that 


--  1  S  -I-  >1  i 

K=^ccr(^-<i»^C:  r  )^ 


1.5  +  A 


CCT 


where 

_3(l-~5 

^ccr)  w  —  ^  ^^CCT  ^^CCT 

"  4(1  +  2  iKccT  +  CcCT  ■ 

Although  this  is  am  implicit  equation  relating  A^  to  A^cr  (and  therefore  must  be 
solved  iteratively),  we  can  see  that  the  equation  is  "gilmost"  explicit  by  noting  that  V^ct 
»  0.2  is  often  a  good  approximation  for  porous  media.  Then  both  Aj  and  Apcr  will  be 
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small  in  magnitude  and,  therefore,  the  last  factor  in  the  right-hand  side  of  the  equation 
will  be  closely  approximated  by  unity.  Thus  is  uniquely  determined  by  and  (j)  in 
the  expected  range  of  the  elastic  parameters  for  porous  materials. 


Step  3  and  4.  SC.  The  moduli  and  of  the  matrix  where  all  inclusions  (of 
concentration  (j))  are  filled  with  the  cement  can  be  found  from  the  following  two 
equations  (Benyman,  1980): 

(1  -  (^)(K^  -  =  0, 

(1  -  ^)(G,  -  +  m  -  G,u,)Qc  =  0, 

where 

p  _  ^Fill  _  ^Fill  ^  Gpiii 

K,+%G„'  K,+%G„' 

Q  _  Gpiii  +  Z  ^  _  Gpjii  -i-  Z  2  _  Gfm(9Kpj,i  -h  8Gf,^) 

G,-hZ’  G,  +  Z’  6(Kp„+2Gp„) 

Next  we  calculate  the  effective  moduli  ( Kp^  and  G^^)  of  the  same  matrix  where  some 
of  the  inclusions  are  filled  with  the  cement,  and  some  are  empty.  If  the  concentration  of 
the  empty  inclusions  is  ,  the  concentration  of  the  cemented  inclusions  is  ^  —  (j)^ 
because  the  volumetric  fraction  of  the  matrix  is  still  1  -  0 .  The  desired  moduli  are  found 
from  the  following  two  equations: 

(1  -  (l»iK^  -  Kpff)P,  +  (0  -  0J(^,  -  Kpff)P,  -  tKp^P,  =  0, 

(1  —  0)(Gj  —  Gpg)Q^  -I-  (0  —  0^)(G^  —  Gpg)Q^  —  ^fGp^Qf)  =  0, 


where 


_  p 


"  YsGe, 


Qs  = 


Gpff  +  Z 

Gs  +  Z' 


Qc  = 


Gsff  ^ 

G,-hZ  ’ 


z  = 


Gpffi^KEff  +  8Gg^) 

^(Kpff+'^Gpff) 
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Steps  3'  and  4'.  DEM.  We  use  the  iterative  scheme  developed  by  Norris  (1985).  We 
choose  a  volumetric  fraction  increment  d(l)  such  that  «  ^ .  The  total  number  of 
iterations  N  is  (Zimmerman,  1991): 


N  =  - 


The  number  of  the  first  steps  is  =  Ncj)^  /  <j) .  During  these  iterations  we 
calculate  the  effective  bulk  and  shear  moduli  as: 

- ^ 

W  M  Ar^  (;) 

GEff^'  ^  =  Gj^^  — ^ - rr^#; 

(0)  _  (0)  _  ^ 


The  last  N  —  iterations  are  performed  using  equations 


=  V  + 


(g,-V^>)(4G,/>  +  3V^>) 

4G,/>  +  3i^, 


Eff  Eff  9/^  o)+8G,/>  +  6(V'^  +  2G^/>)G,/G^/>  ^■ 
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APPENDIX  2:  COMPARISON  OF  VARIOUS  EFFECTIVE  MEDIUM  SOLUTIONS 


Below,  we  compare  the  results  of  various  effective  medium  theories  (see  text  for  the 
meaning  of  abbreviations)  for  two-phase  composites  made  of  (a)  glass  £ind  epoxy,  eind  (b) 
quartz  and  ice.  The  verticcd  lines  in  the  graphs  correspond  to  the  volumetric  fraction  of 
0.36  which  is  approximately  the  fraction  of  cement  in  the  case  where  the  spherical 
inclusions  touch  each  other.  Only  the  Lower  HS  limit  comes  close  to  the  experimental 
results  and  our  new  solution  (see  Figure  4). 


Figure  7.  Compresslonal  and  shear  moduli  of  two-phase  composites  as  given  by  various 
effective  medium  theories. 
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2.2.  DIAGNOSING  SANDSTONES:  STRENGTH  AND  PERMEABILITY  FROM  POROSITY 

AND  VELOCITY 

In  collaboration  with  Ivar  Brevik 

INTRODUCTION 

Non-uniqueness  in  relating  veloeity  to  porosity  in  core  and  well-Iog  data  complicates 
interpretation  of  sonie  and  seismic  measurements.  One  reason  for  this  non-uniqueness 
in  sandstones  is  elay  (e.g.,  Han,  1986).  Another  reason  is  textural  variability  among 
samples.  Dvorkin  and  Nur  (1996)  examine  two  relatively  elay-free  sandstone  groups  in 
the  same  porosity  range,  but  whose  veloeities  significantly  differed  (Figure  la).  By 
eomparlng  the  data  with  effective-medium  theories  they  Interpret  this  velocity  difference 
as  resulting  from  the  difference  in  the  position  of  diagenetle  eement.  The  explanation  is 
that  in  the  "fast"  (Oseberg)  rocks  (contact)  cement  is  located  predominantly  at  the  grain 
eontacts  whereas  in  the  "slow"  (Troll)  rocks  (non-contact)  cement  is  loeated 
predominantly  away  from  these  contacts. 

Coincidentally,  the  permeability  of  the  Troll  roeks  is  smaller  than  that  of  the  Oseberg 
rocks  (Figure  lb).  This  faet  allows  us  to  assume  that  the  position  of  the  diagenetie 
eement  affeets  not  only  veloeity  but  also  permeability.  The  assumption  is  supported  by 
numerical  simulations  of  Bosl  et  al.  (1998).  This  effect  has  a  simple  physieal 
explanation:  the  non-eontact  cement  acts  to  increase  the  specific  surface  area  (Figure 
lb)  and  thus  decrease  permeability. 

Consider  now  a  dataset  where  a  clear  trend  between  porosity  and  permeability  is 
absent.  Our  working  hypothesis  is  that  this  may  be  due  to  the  variability  in  the  loeation 
of  diagenetie  cement.  We  use  this  hypothesis  in  a  case  study  where  permeability  and 
porosity  data  are  available  from  cores  in  a  well  and  velocity  Is  available  from  a  sonie  log. 
A  relationship  between  permeability  and  porosity  is  absent.  The  goed  is  to  find  a  textural 
property  of  high-porosity  samdstones  with  whieh  permeability  eorrelates  well. 
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To  do  so,  we  use  log-derived  velocity  and  porosity  together  with  rock  physics  theory 
to  subdivide  the  diagenetic  cement  Into  the  contact  and  non-contact  parts.  Then  by 
relating  permeabiUty  to  the  volumetric  fraction  of  the  non-contract  cement  in  the  rock 
we  obtain  a  usable  trend.  This  non-contact  cement  fraction  is  the  desired  textural 


property. 


Figure  1.  a.  Compressional  modulus  (the  product  of  bulk  density  Eind  P-wave  velocity 
squared)  versus  porosity  for  the  Oseberg  and  Troll  samples.  The  data  displayed  are  for 
room-diy  rocks  at  30  MPa  effective  pressure  (Strandenes,  1991;  and  Blangy,  1992).  The 
upper  curve  is  from  the  contact  cement  theory  and  the  bottom  curve  is  from  the  non-contact 
cement  theory  (Dvorkin  and  Nur,  1996).  b.  Permeability  versus  porosity  for  the  same 
datasets  (Strandenes,  1995).  The  cartoons  schematically  show  the  location  of  cement 
among  grains  (contact  for  Oseberg  and  non-contact  for  Troll). 


By  quantifying  the  amount  of  contact  and  non-contact  cement  we  diagnose 
sandstone's  texture  from  well  log  data,  similar  to  Dvorkin  and  Nur  (1998).  We  show  that 
such  diagnostic  is  important  not  only  for  obtaining  a  high-correlation  permeability  trend 
but  also  for  assessing  the  strength  of  the  rock. 


NORTH  SEA  SLEIPNER  FIELD,  WELL  15/9-16 

A  vertical  well,  15/9-16,  penetrates  the  North  Sea  Sleipner  gas/condensate  reservoir 
comprised  of  Paleocene  turbiditic  sand.  Porosity  emd  permeability  are  available  from 
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about  60  plugs  that  evenly  cover  the  Interval  from  2380  to  2460  m.  The  verUcal  and 
horizontal  pemteabUiUes  are  practically  Identical.  The  latter  is  used  In  this  study.  The 
sandstone  Is  veiy  well  sorted  with  the  average  sorting  coefficient  of  1.7;  grain  size  vanes 
between  0.15  and  0.25  mm.  The  grains  are  predominantly  quartz  (average  80%)  with  the 
rest  being  feldspar  (average  14%).  mica  (average  2.3%).  and  clay.  mosUy  chlorite,  (average 
2.2%).  Traces  of  calclte  and  other  minerals  are  also  present.  The  contact  cement  in 
these  rocks  is  quartz  (Nadeau,  1998). 

The  upper  part  of  the  well  is  saturated  with  gas,  with  the  gas-water  contact  at  2430 
m.  The  bulk  moduli  and  densities  of  the  formation  water  and  gas  are  2.75  GPa  and  1.02 
g/cm^,  and  0.07  GPa  and  0.27  g/cm^,  respectively  (calculated  using  Batzle  and  Wang, 

1992). 

The  interval  under  investigation  can  be  subdivided  into  a  high-resistivity  zone  (HR^ 
overlaying  a  low-resistivity  zone  (LRZ)  with  the  transition  at  about  2410  m  (Figure  2a). 
Nadeau  (1998)  shows  that  there  is  a  diagenetic  change  associated  with  this  transition. 
HRZ  has  a  restricted  distribution  of  diagenetic  chlorite  and  up  to  5%  quartz  cement. 
LRZ  bag  a  slightly  larger  content  of  chlorite  and  a  smaller  degree  of  cementation. 

Porosity  in  the  entire  interval  is  high.  We  calculate  it  from  bulk  density 
(Schlumberger,  1989).  It's  values  do  not  differ  much  from  those  measured  on  cores 
except  for  a  few  points  in  LRZ  (Figure  2b).  The  final  results  of  this  study  practicaUy  do 
not  depend  on  what  porosity  (log-derived  or  core)  we  use.  For  this  reason  we  relate  all 
parameters  to  the  log-derived  porosity. 

The  available  log  data  contain  only  compressional-wave  velocity.  We  calculate  the 
dry-frame  compressional  modulus  (the  product  of  bulk  density  and  P-wave  velocity 
squared)  by  the  V^-only  fluid  substitution  method  of  Mavko  et  al.  (1995).  The  result  is 
well  matched  by  the  values  measured  on  several  selected  room-dry  plugs  (Figure  2d).  The 
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30  MPa  effective  pressure  for  these  datapoints  essentially  equals  the  reservoir  effective 
pressure  (about  29  MPa). 


1  10  100  1000 
C  Permeability  (mD) 


j  10  20  30 

U  Comp.  Modulus  (GPa) 


Figure  2.  Various  parameters  versus  depth  in  well  15/9“16.  a.  Far  resistivity.  Gray  curve 
is  for  LRZ.  b.  Log-derived  (gray  curve)  and  core  porosity.  Open  S3mibols  are  for  LRZ.  c. 
Permeability.  Open  symbols  are  for  LRZ.  d.  Dry-rock  (gray  curve)  and  directly  measured 
(black  curve)  compressional  modulus.  Symbols  are  from  dry-rock  lab  measurements  at  30 
MPa. 


ANALYSIS  OF  DATA 

There  is  a  very  weak  and  practically  not  useful  correlation  between  permeability  and 
porosity  (either  log-derived  or  core)  in  well  15/9-16  (Figure  3).  This  fact  can  also  be 
observed  in  Figures  2b  aind  2c  where  porosity  and  permeability  are  plotted  versus  depth. 
The  lower-porosity  sandstones  in  HRZ  have  permeability  larger  than  that  of  the  higher- 
porosity  sandstones  in  LRZ.  Notice  also  that  the  LRZ  sandstones  are  softer  than  the 
HRZ  sandstones  (Figure  2d). 

The  modulus-porosity  and  permeability-porosity  trends  for  HRZ  and  LRZ  cire  given  in 
Figure  4.  In  the  modulus-porosity  plane  (Figure  4a)  the  HRZ  trend  parallels  that  of  the 
contact-cemented  Oseberg  rocks  and  the  contact  cement  theoretical  trajectoiy  (the  latter 
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calculated  for  quartz  grains  with  quartz  cement).  The  HRZ  rocks  have  a  tight  modulus- 
porosity  trend.  Remarkably,  in  the  permeability-porosity  plane  (Figure  4b)  the  HRZ 
sandstones  plot  on  top  of  the  Oseberg  data  and  also  exhibit  a  noticeable  permeability- 
porosity  trend. 


Figure  3.  Horizontal  permeability  versus  log-derived  (a)  and  core  (b)  porosity.  Gray  lines 
show  best  linear  fits.  Correlation  coefficients  are  given  in  the  graphs. 


Figure  4.  a.  Dry-frame  compressional  modulus  versus  porosity  for  the  Oseberg  and  Troll 
samples  at  30  MPa,  and  HRZ  and  LRZ.  The  HRZ  and  LRZ  data  are  selected  at  the  depths  of 
permeability  datapoints.  The  upper  curve  is  from  the  contact  cement  theory  and  the  bottom 
curve  is  from  the  non-contact  cement  theory,  b.  Permeability  versus  porosity  for  the  same 
datasets.  The  open  triangles  are  for  Oseberg  and  Troll.  The  filled  circles  are  for  HRZ  and 
the  open  circles  are  for  LRZ.  The  Oseberg  and  Troll  data  are  plotted  versus  core  porosity 
whereas  the  HRZ  and  LRZ  data  are  plotted  versus  log-derived  porosity. 
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The  LRZ  sandstones  on  the  other  hand  do  not  show  a  modulus-porosity  trend. 
These  datapoints  fill  the  space  between  the  contact  and  non-contact  cement  theoretical 
trajectories  (the  latter  calculated  for  quartz  rock  at  29  MPa  effective  pressure).  Similarly, 
in  the  permeability-porosity  plane  these  rocks  fill  the  space  between  the  Troll  and 
Oseberg  datapoints,  and  a  permeability-porosity  trend  is  absent. 

As  we  mentioned  before,  both  HRZ  and  LRZ  rocks  are  quartz-cemented  and  have  very 
close  mineralogy.  Based  on  these  facts,  we  assume  that  both  the  modulus-porosity  and 
permeability-porosity  non-uniqueness  (Figure  4)  is  due  to  the  varying  amounts  of 
contact  and  non-contact  cement. 

DIAGNOSING  ROCK  FOR  NON-CONTACT  CEMENT 

In  order  to  calculate  the  amount  of  the  non-contact  cement  from  sonic  cind  porosity 
we  use  a  model  where  a  high-porosity  sandstone  has  an  idealized  texture.  Its  basic 
framework  is  a  random  dense  pack  of  identical  spherical  grains  at  some  critical  porosity 
(j)^  which  may  vary  between  0.36  and  0.4  (Nur  et  al.,  1998).  Every  grain  is  identiccdly  and 
evenly  enveloped  by  a  layer  of  contact  cement,  and  the  rest  of  the  solid  phase  (additional 
to  the  greiins  and  contact  cement)  forms  non-contact  cement  deposited  in  the  pore  space 
away  from  grain  contacts  (Figure  5a). 

One  is  unlikely  to  encounter  such  an  idealized  picture  in  a  thin  section.  However, 
this  is  a  way  to  build  effective  medium  models,  some  of  which  have  been  effectively  used 
to  model  granular  rocks  (e.g.,  Dvorkin  and  Nur,  1996;  Moos  et  al.,  1997). 

Let  us  now  consider  a  datapoint  in  the  modulus-porosity  plane  that  lies  below  a 
theoretical  contact  cement  trajectory  (Figure  5b).  We  assume  that  the  non-contact 
cement  does  not  contribute  to  the  stiffness  of  this  rock.  Therefore,  its  elastic  modulus  is 
identical  to  that  of  a  higher  porosity  (0^^^)  rock  of  the  same  texture  but  without  the 
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non-contact  cement.  The  corresponding  datapolnt  is  the  horizontal  projection  (in  the 
modulus-porosity  plane)  of  the  original  one  onto  a  contaet-cement  trajectory  (Figure  5b). 


Non- 


Contact  Contact  Cement 


a 


Figure  5.  a.  Idealized  picture  of  granular  rock  with  contact  and  non-contact  cement,  b. 
Calculating  the  amount  of  non-contact  cement.  Gray  s3mibols  show  some  scattered 
datapoints.  c.  Volumetric  fraction  of  contact  (solid  line)  and  non-contact  (dotted  line) 
cement  versus  depth  in  well  15/9-16.  Gray  vertical  bar  shows  the  extent  of  HRZ. 


Volume  balance  gives  the  following  relations  between  the  porosity  of  the  datapoint 
(,(p),  its  projection  on  the  contact-cement  trajectory  critical  porosity  (0^),  and  the 

volume  fractions  of  the  contact  (/^^)  and  non-contact  (f„cc)  cement  in  rock: 

fee  =  -  <Peem  ^  fnee  =  ^  “  /cc*  ^ 

The  contact-cement  trajectory  can  be  plotted  using  equations  in  Dvorkin  and  Nur 
(1996):  it  depends  on  the  elastic  moduli  of  grains  and  cement,  and  on  the  chosen  critical 
porosity  value.  An  additional  input  parameter  is  the  average  number  of  contacts  per 
grain  in  the  original  sphere  pack  (n).  It  may  vary  between  9  and  8. 

In  the  case  under  examination  we  choose  =  0.38  (average  between  0.36  and  0.4) 
and  n  -  8.5.  Because  the  grains  and  contact  cement  are  predominantly  quartz,  we  plot 
the  contact-cement  trajectory  for  quartz-cemented  quartz  grains  (Figure  5b).  In 
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calculating  this  trajectory,  we  use  38  GPa  and  44  GPa  for  the  bulk  and  sheeir  moduli  of 
quartz,  respectively  (Carmichael,  1990). 

This  trajectory  can  be  statisticcdly  fitted  (with  correlation  coefficient  about  1)  by  the 
equation 

1.037-^0.38 = -0.0013 +  0.0134M  + 4.3- 10"^  (2) 

where  M  is  the  compresslonal  modulus  of  the  dry  contact-cemented  rock.  The  contact- 
cement  trajectory  for  emy  set  of  input  parameters  can  be  obtained  from  Dvorkln  and  Nur 
(1996).  Then  it  can  be  statistically  fitted  by  an  equation  similar  to  Equation  (2). 

Now  we  can  combine  Equations  (1)  and  (2)  to  arrive  at  a  formula  that  relates  the 
volumetric  fractions  of  the  contact  and  non-contact  cement  to  the  dry-frame 
compresslonal  modulus  and  porosity: 

=  0.93(-0.0013-h  0.0134M -1-4.3  •  10-' =  0.38  -  0  (3) 

These  fractions,  as  calculated  for  well  5/9-16  (using  the  dry-frame  compresslonal 
modulus  from  fluid  substitution)  are  given  in  Figure  5c.  The  contact  cement  dominates 
in  HRZ  whereas  the  non-contact  cement  is  prevalent  in  LRZ. 

This  diagnostic  can  be  immediately  used  to  assess  the  strength  of  the  rock:  clearly 
the  larger  the  amount  of  contact  cement  the  stronger  the  rock  (at  the  same  porosity). 
This  effect  could  be  clearly  seen  in  the  Troll  and  Oseberg  example  (Figure  1).  The  TroU 
samples  that  do  not  have  contact  cement  are  friable  sands  (Blangy,  1992),  whereas  the 
Oseberg  samples  show  significant  structured  integrity  (Stremdenes,  1991).  In  the  case 
under  examination,  our  diagnostic  is  also  consistent  with  the  rock's  strength:  Nadeau 
(1998)  states  that  quartz  cementation  progressed  more  readily  in  HRZ  emd  is  associated 
with  intervals  less  prone  to  semd  production. 
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PERMEABILITY  TREND 


In  Figure  6a  we  plot  the  logarithm  of  permeability  (k)  versus  the  volumetric  fraction 
of  the  non-contact  cement.  A  linear  trend  is  evident  (as  opposed  to  the  absence  of  such 
in  Figure  3).  The  linear-fit  equation  for  this  trend  is 

Log,o^  =  3.3  - 19.46/„,,;  R  =  0.85.  (4) 

The  correlation  slightly  improves  if  the  permeability  is  normalized  by  the  grain  size 
(d)  squared  (Figure  6b): 

Log,o(it  /  J')  =  4.8  -  20.47/^,;  R  =  0.86.  (5) 

Such  normalization  is  often  used  to  improve  permeability  trends  because  permeability 
strongly  depends  on  the  grain  size  (e.g.,  Bourbie  et  aL,  1987).  In  our  case  this 
improvement  is  small  due  to  a  relatively  uniform  grain  size  distribution  in  the  interval 
(between  0.15  aind  0.25  mm). 


Figure  6.  Permeability  (a)  and  permeability  normalized  by  grain  size  squared  (b)  versus  the 
volumetric  fraction  of  non-contact  cement.  Gray  lines  show  best  linear  fits.  Correlation 
coefficients  are  given  in  the  graphs. 

CONCLUSION  AND  RECOMMENDATIONS 

The  high-correlation  trends  expressed  by  Equations  (4)  and  (5)  can  be  used  to  predict 
permeability  from  sonic  and  porosity  data.  We  obtained  these  trends  by  using  the 
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amount  of  non-contact  cement  to  quantify  the  variability  of  roek  texture.  This  eonfirms 
our  initial  hypothesis  that  textural  variability  in  rock,  specifically,  the  amount  of  the 
non-contact  cement,  can  be  responsible  for  non-uniqueness  in  permeability-porosity 
trends  sueh  as  shown  in  Figure  3.  To  ealeulate  the  amount  of  the  non-contaet  cement 
is  a  way  of  obtaining  a  meaningful  eorrelation  for  permeability. 

The  method  of  ealculating  the  amount  of  contact  and  non-contact  cement  presented 
here  is  in  fact  a  method  of  diagnosing  the  texture  of  high-porosity  sandstone  from  well- 
log  data.  Such  diagnostic  is  important  not  only  for  obtaining  a  usable  correlation  for 
permeability  but  also  for  assessing  the  strength  of  rock  and  its  susceptibility  to  Scinding. 

Equation  (2)  is  not  universal.  In  every  separate  ease  study  it  has  to  be  re-derived 
from  the  appropriate  contact-eement  trajeetory  where  the  elastic  moduli  of  the  grains 
and  cement  are  selected  aeeording  to  mineralogy  or  within  reasonable  ranges.  The  "free" 
parameters  (j)^  and  n  can  be  varied  within  reasonable  (narrow)  ranges  until  the  best 
correlation,  such  as  Equation  (4),  is  obtained. 
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2.3.  DIAGNOSING  SANDS  USING  SONIC  AND  SEISMIC 
In  collaboration  with  Per  Avseth,  Gaiy  Mavko,  and  Johannes  I^kkje 

ABSTRACT 

At  high  porosity,  velocity  in  reservoir  rocks  strongly  depends  on  the  position  of  the 
intergranular  material.  Velocity  is  high  if  the  original  grains  are  cemented  at  their 
contacts.  It  is  low  if  the  pore-fllling  material  is  placed  away  from  the  contacts.  In  the 
latter  case  we  have  truly  unconsolidated  sediments.  In  the  former  case  we  have  high- 
porosity  cemented  rocks.  Separating  these  two  rock  types  is  important  for  hydrocarbon 
identification.  Due  to  the  difference  in  the  rock  frame  stiffness  between  the 
unconsolidated  and  high-porosity  cemented  rocks,  seismic  signatures  of  the  former  filled 
with  water  can  be  very  close  to  those  of  the  latter  filled  with  hydrocarbons.  This  may 
complicate  direct  hydrocarbon  detection.  We  separate  the  two  rock  types  by  diagnosing 
sand  using  rock  physics  theory.  We  conduct  such  diagnostic  on  well  log  data  from  two 
wells  that  penetrate  the  Heimdal  formation  (North  Sea).  We  show  that  the  Heimdal 
formation  reservoir  is  composed  of  both  unconsolidated  and  cemented  high-porosity 
sands.  The  initial  quartz  cementation  present  in  the  latter  is  clearly  seen  in  the 
cathode-luminescent  SEM  images.  These  images,  combined  with  point  XRD  analysis, 
confirm  our  diagnostic  that  the  high-velocity  high-porosity  sands  in  Heimdal  contain 
quartz  grains  surrounded  by  quartz-cement  rims.  We  find  that  the  two  different  types  of 
sand  which  are  capped  by  similar  low-impedance  shales  produce  drastically  different 
AVO  signatures.  The  oil-filled  high-porosity  cemented  sand  shows  a  relatively  strong 
zero-offset  reflectivity  which  becomes  less  positive  wdth  increasing  offset,  while  the  oil- 
filled  uncemented  sand  shows  a  negative  zero-offset  reflectivity  wdth  Increasingly 
negative  far-offset  response.  These  results  show  that  (1)  rock  diagnostic  can  be 
conducted  not  only  on  the  log  scale  but  edso  on  the  seismic  scale;  and  (2)  taking  into 
account  the  nature  of  the  rock  improves  our  ability  to  identify  pore  fluid  from  seismic. 
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INTRODUCTION 


Quartz  cementation  of  sands  greatly  affects  porosity,  permeability,  and  seismic 
properties.  Sandstones  in  continuously  subsiding  sedimenteiry  basins,  such  as  in  the 
North  Sea  and  the  Gulf  Coast,  tend  to  have  poorly  developed  quartz  cement  down  to  a 
depth  of  2.5  -  3.0  km  (Bjorlykke  and  Egeberg,  1993).  Hence,  Tertiary  sands  in  the  North 
Sea  are  usually  reported  to  be  poorly  consolidated  with  no  (or  insignificant  quantities  of) 
quartz  cement.  “Insignificant”  is  related  to  volume  —  small  amounts  of  quartz  cement 
do  not  significantly  affect  porosity.  However,  only  small  amounts  of  cement  at  grain 
contacts  are  needed  to  considerably  stiffen  the  frame  of  a  rock  (Dvorkin  and  Nur,  1996) 
and  strongly  increase  velocity.  We  apply  the  contact  cement  concept  to  study  two  clean 
sandstone  intervals,  both  representing  the  Palaeocene  age  Heimdal  Formation  in  the 
North  Sea.  Both  intervals  are  oil-filled  reservoir  sands  of  commercial  Interest.  We 
diagnose  the  rock  using  well  log  measurements  and  rock  physics  theory.  We  assume 
that  If  in  the  velocity-porosity  plane  a  datapoint  falls  close  to  a  theoretical  line,  the 
internal  structure  of  the  rock  is  similar  to  the  idealized  structure  used  in  the  model.  We 
find  from  such  diagnostic  that  one  interval  is  composed  of  unconsolidated  sand,  while 
the  other  interval  is  composed  of  cemented  high-porosity  sand.  Thin-sectlon  and  SEM 
images  confirm  this  diagnostic.  By  studying  the  seismic  signatures  of  these  two  different 
types  of  clean  sands  we  upscale  the  log-based  diagnostic  to  the  seismic  scale. 

ROCK  DIAGNOSTIC  AND  CONFIRMATION 

We  examine  two  wells  —  Well  #1  and  Well  #2.  Sonic  velocity  and  gamma-ray  are 
plotted  versus  depth  for  both  wells  in  Figures  la  to  Id.  Vp  is  plotted  versus  porosity  in 
Figures  Id  and  If.  Notice  that  in  Well  #2  a  thick  sand  interval  (gray  bar  in  Figure  Ic)  is 
marked  by  extremely  low  and  constant  gamma-ray  readings.  This  sand  layer  is 
surrounded  by  high-gamma-ray  shale  packages.  In  Figure  If,  these  two  lithologies  fall 
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into  two  distinctive  velocity-porosity  patterns.  In  Well  #1,  unlike  in  Well  #2,  we  observe 
a  gradual  variation  of  clay  content  between  very  cleatn  sand  and  shale.  Only  a  relatively 
thin  (10  m)  sand  interveil  (gray  bar  in  Figure  la)  is  identified  as  a  practically  clay-free 
reservoir  sand.  Because  of  the  gradual  variation  of  clay  content  in  this  well,  we  do  not 
observe  (Figure  le)  velocity-porosity  patterns  as  distinctive  as  in  Well  #1.  These  two 
clean  sand  intervals  (in  both  wells)  represent  the  same  stratigraphic  level,  although 
located  in  different  oil  fields.  They  are  shown  by  bold  black  symbols  in  Figures  le  and  If. 


2  3  40  80  120  2  3  4  40  80  120 


Figure  1.  P-wave  velocity  and  gamma-ray  versus  depth  (a-d);  and  P-wave  velocity  versus 
porosity  (e  and  f)  for  both  wells. 

To  directly  confirm  this  diagnostic,  consider  the  thin  sections  of  two  samples  from 
both  intervals  (Figure  3).  The  samples  have  approximately  same  porosity.  They  are 
predominantly  queirtz.  No  contact  cementation  is  appairent  in  both  images.  The  left 
image  (Well  #1)  shows  some  organic  coating  around  quartz  grains.  Consider  now  two 
SEM  Images  of  a  sample  from  Well  #2  (Figure  4).  The  left-hand  image  is  in  back-scatter 
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light  and  the  right-hand  one  is  in  eathode-luminescent  light.  Notiee  the  V-shaped  grain 
in  the  middle.  No  contact  cement  rim  is  apparent  around  this  grain  in  back-scatter 
light.  Cathode-luminescent  light  reveals  a  contact-cement  rim  around  this  grain.  The 
point  XRD  analysis  shows  that  both  the  grain  and  cement  rim  are  pure  quartz.  This 
confirms  our  diagnostic  that  the  Well  #2  sand  interval  is  contact-cemented.  The 
hexagonal  crystal  shapes  in  the  upper  left  comer  cdso  indicate  diagenetic  cementation. 
No  cement  rims  or  hexagonal  crystal  shapes  have  been  found  around  grains  in  the  sand 
Interval  from  Well  #1.  Another  direct  proof  of  our  diagnostic  was  that  cores  extracted 
from  Well  #  1  appeared  as  piles  of  loose  sand,  whereas  those  from  Well  #2  supported 
external  stress. 


SEISMIC  RESPONSE 

To  understand  how  the  type  of  sand  (unconsolidated  versus  cemented)  affects  the 
seismic  response,  we  analyze  CDP  gathers  at  the  well  locations.  Figure  5a  shows  the 
real  CDP  gather  at  Well  #1  where  the  picked  horizon  is  at  the  top  of  the  Heimdal 
formation.  Figure  5b  gives  a  S3nithetic  CDP  gather  for  this  well  produced  by  using  a  30 
Hz  zero-phase  Ricker  wavelet  and  a  log-derived  reflectivity  series.  Both  the  reed  and 
synthetic  gathers  show  reflectivity  increasingly  negative  with  increasing  offset  at  the 
picked  horizon.  This  reflectivity  is  plotted  versus  offset  (angle),  together  with  the 
theoretical  Zoepprltz  line,  in  Figure  5e.  Contrary  to  Well  #1,  the  top  of  the  Heimdal 
formation  in  Well  #2  (which  is  capped  by  similar  shales)  produces  a  strong  positive 
reflector  with  reflectivity  decreasing  with  increasing  offset  (Figures  5c  and  5d).  For  this 
well,  the  reflectivity  is  plotted  versus  offset  (angle),  together  with  the  theoretical  Zoepprltz 
line,  in  Figure  5f.  The  S5mLthetic  response  is  very  close  to  the  real  data  in  both  wells 
which  means  that  we  can  rely  on  well-log-based  rock  diagnostic  to  predict  seismic 
response. 
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Figure  2.  P-wave  velocity  versus  porosity  for  sand  intervals  in  both  wells.  Theoretical 
lines  serve  to  diagnose  the  rocks. 


Figure  3.  Thin  sections  of  two  selected  samples  from  Well  #1  (left)  and  #2  (right). 


Figure  4.  SEM  images  of  a  Well  #2  sample  in  back  scatter  light  (left)  and  cathode- 
luminescent  light  (right). 


This  offset-dependent  reflectivity  analysis  shows  that  clean  sands  of  the  same 
formation,  similar  porosity,  and  with  comparable  oil  saturation  produce  drastically 
different  seismic  response  depending  on  whether  they  are  truly  unconsolidated  or  have 
initial  quartz  cementation.  Therefore,  we  can  use  both  normal-incidence  and  offset- 
dependent  reflectivity  to  diagnose  rock  and  characterize  a  reservoir  from  seismic.  Such 
rock  diagnostic  may  be  of  great  importance  because  if  high-porosity  contact-cemented 
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sands  are  not  separated  from  truly  unconsolidated  sands,  one  may  misinterpret  a 
change  in  seismic  signatures  caused  by  this  petrographic  effect  as  a  pore-fluid  effect. 


Figure  5.  Top.  Real  (a  and  e)  and  synthetic  (b  and  d)  CDP  gathers.  In  synthetic  gathers, 
the  AVO  effect  was  modeled  only  at  the  target  zones.  Bottom.  Real  reflectivity  versus  offset 
and  angle  (symbols)  and  theoretical  Zoepprltz  lines. 


CONCLUSION 

Rock  diagnostic  is  important  for  correctly  characterizing  prospective  reservoirs.  Such 
diagnostic  is  based  on  rock  physics  theory  and  can  be  accomplished  using  well  log  data. 
The  diagnostic  features  observed  in  well  log  data  can  be  translated  into  distinctive 
seismic  signatures.  Therefore,  seismic  data  can  edso  be  used  for  rock  diagnostic  given 
that  the  stratigraphic  unit  is  resolvable  at  the  seismic  scale.  In  this  paper  we  applied  the 
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diagnostic  concept  to  the  Heimdal  formation,  eind  were  able  to  discriminate  high- 
porosity  cemented  from  unconsolidated  sands  both  from  well  logs  eind  seismic. 
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